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Abstract. We study a simple model of a bouncing ball that takes explicitely 
into account the elastic deformability of the body and the energy dissipation 
due to internal friction. We show that this model is not subject to the problem 
of inelastic collapse, that is, it does not allow an infinite number of impacts in 
a finite time. We compute asymptotic expressions for the time of flight and 
for the impact velocity. We also prove that contacts with zero velocity of the 
lower end of the ball are possible, but non-generic. Finally, we compare our 
findings with other models and laboratory experiments. 



1. Introduction 

In this paper we study how a ball bouncing against an horizontal, rigid plane 
comes to rest. We are motivated by the fact that bouncing objects are the basic 
building blocks of granular fluids. Our view is that an acceptable mathematical 
description of granular fluids is impossible without taking into account internal 
vibrations of the bouncing objects. To explore the viability of this program we 
formulate a simple, one-dimensional model that explicitly accounts for the defor- 
mation of the ball, and we prove that, while encompassing the desirable properties 
of current models, which neglect internal vibrations, it does not incur their patholo- 
gies. 

The simplest and most widely used model of a bouncing ball (or grains of a 
granular fluid) assumes that the ball is a rigid body, and that an impact with 
the floor is an instantaneous event, which reverses the vertical component of the 
speed of the ball. In order to model energy dissipation caused by an impact, it 
is customary to introduce a positive coefficient of restitution r < 1, so that the 
vertical speed u a immediately after an impact is related to the vertical speed u 
immediately before the impact by the simple relationship 

(1) u a = -ru b . 

This model performs well in cases where the ball does not experience too many 
impacts in the unit of time. It has been used, for example, as an ingredient in 
the description of sport balls [I], and to study the dynamics of a bead on a vi- 
brating plate (for its deep mathematical facetes the latter has become a classical 
problem, see [2] sec. 2.4). The assumptions behind this model (namely: neglect of 
deformability, instantaneous impacts and energy losses described by a restitution 
coefficient) are the basic building blocks of current theories of granular gases [5] . 

The most apparent drawback in this approach is that it cannot limit a priori 
the number of impacts in the unit of time. In fact, granular systems described 
by J!}, for a large class of parameter choices and initial conditions, are subject to 
the phenomenon of inelastic collapse, where clusters of particles are subject to an 
infinite number of collisions in a finite time [I]. The simplest example of inelastic 
collapse is given by a single ball, subject to a constant gravity force, bouncing 
repeatedly off the floor. Neglecting the interaction with the air, the vertical speed 
of the ball immediately after the (n + 1)— th impact is linked to that at the n— th 
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impact by 

(2) u n+ i = ru n . 

The duration of the n— th flight is r„ = 2u n /g, where g is the acceleration of gravity. 
The sum of the times of flight is easily computed, and gives a geometric series that 
converges to 

_ 2u 
5(1 - r) 

For times larger than the model is meaningless. 

This pathology is not caused by the one-dimensional nature of the model sketched 
here. Even taking into account rotational degrees of freedom and exchanges of 
angular momentum at the impacts due to surface friction, the inelastic collapse is 
still a common outcome [5]. 

To avoid the problem it is customary to assume that the restitution coefficient 
is an increasing function of the impact velocity, and that r(u) — > 1 if u — > 0. A 
popular choice is 




where U is a constant that depends on the material and geometry of the ball. 
This expression has some theoretical support [6], and it appears to fit the data 
for non- repeated impacts [J]. In the absence of gravity, inelastic collapse is ruled 
out rigorously for systems of three balls [3], and there is numerical evidence that 
the same result holds for an arbitrary number of particles (see [5] chap. 26 and 
references therein). However, using (J3j) in (2]) , yields again a converging sum of 
the times of flight for the bouncing ball problem. This is proved rigorously in the 
appendix (Remark IQ.2|) . 

The common wisdom is that is the time when the bouncing ball comes to 
rest and maintains permanent contact with the floor (see, for example [3 J chap. 3). 
This point of view faces serious difficulties if one desires to model, for example, 
the settling-down of many beads poured into a box: it is not clear how a bead 
at rest on the floor should behave when hit by a moving one, because there is no 
obvious way to extend (p} to situations where three or more bodies are in contact 
simultaneously. 

More generally, the assumption (TJ rests on hypoteses that become invalid as 
the frequency of the impacts diverges: collisions with the floor are not truly in- 
stantaneous, and treating the ball as indeformable is highly questionable when the 
frequency of impacts is close to the resonant frequencies of the bouncing ball (e.g. 
[9] sec. 2.2). 

Of course, a model avoiding inelastic collapse must at the same time account for 
the obvious observation that bouncing balls do come to rest after a finite time. A 
first step in this direction is the time-of-contact model [10], which prescribes r < 1 
if the time-of-flight of the ball is greater than a constant T c , and r = 1 otherwise: 
conservative impacts are interpreted as internal vibrations of the ball, which (from 
a macroscopic point of view) maintains contact with the floor. A variation on this 
theme leads to the notion of a stochastic restitution coefficient [TT] . 

In this paper we study a one-dimensional model of a bouncing ball simple enough 
to allow for rigorous mathematical analysis, but including all the elements that we 
believe are important for further developments in the description of granular fluids. 
We explicitly take into account the deformability of the body, and at the same time 
we give up the notion of restitution coefficient, at least as a primitive concept. We 
prove that our model is free from pathologies analogous to the inelastic collapse, and 
that it is well defined at all times. We still assume impacts to be instantaneous, but 
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Figure 1 . A bouncing ball of mass 2m and diameter L is idealized 
as two points of mass m connected by a massless, dissipative spring. 

only from a microscopic point of view: a persistent contact with the floor is seen as 
a rapid sequence of instantaneous impacts. This property makes the dynamics only 
piecewise-smooth. There is a growing literature on piecewise smooth dynamical 
systems, with engineering and plasma physics applications (for a review see |12|). 
However, the focus is generally on periodically forced systems, rather than on the 
pathways chosen by an unforced system to reach the asymptotic equilibrium. 

The rest of the paper is organized as follows: in section [2] we describe the model 
and we state the main results: in particular the absence of inelastic collapse and 
how the system approaches the rest state; the theorems are proved in sections [3] 
through [6j in section [7| we show that the notion of restitution coefficient is naturally 
recovered, as a consequence of the dynamics, when the time-of-flight is large with 
respect to the characteristic damping time of the internal vibrations; numerical 
simulations are compared with the laboratory experiments of [13]; finally, section [8] 
contains a summary of the results and some forward-looking remarks. 

2. Equations of Motion 

An basic model of a deformable bouncing ball is shown in figure ([I]) : two point 
masses are connected by a massless dissipative spring. This idealized ball, when it 
is not in contact with the floor, is ruled by the following equations of motion 

mx = —gm + k(y — x — L) + v (y — x) 
my = —gm — k(y — x — L) — v{y — x) 

where m is the mass of the material points at x(t) and y(t), L is the length at rest 
of the spring, k is its elastic constant, v is a damping coefficient, and —g is the 
acceleration of gravity. All constants are positive. 

We indicate by t n the time when an impact occurs, i.e. x(t n ) = 0. We define 
the time of flight 

(5) T n — t n -\-\ t n 

and we assume to = 0. Impacts are modeled as an instantaneous elastic interaction 
obeying to the rule 

(6) i(t+) = -z(t-). 

where the notation f{a il ) means lim t ^ a ± fit). The boundary does not exert any 
force directly on the mass in y, which is affected by boundary hits only through 
the resulting deformation of the spring. In other words, at the impact times t n the 
upper point mass obeys the rule 

(7) y(t+) = y(t-). 
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The special (t n ) — may lead to a continuous contact of the lower point 

mass with the floor. This is discussed in detail in the next section. 

We re- write the equations (0]) using L as the scale of lengths, and \Jmj (2k) as 
the scale of times. The resulting dimensionless equations are 



(8) 



x + fix + \x = -7 - \ + \y + fiy 
y + fiy + = — 7 + | + + fix 



where 7 = gm/(2Lk) and /1 = v j\/2km. 

The positions x(t) and y(t) are defined for t g [0, too), where too is the time 
of inelastic collapse, should it occur, or 00 otherwise. For t e U n (t n ,t n+1 ), the 
equations ((8j) guarantee the existence of the derivatives x^ and y^ k > for all k. For 
brevity, it is convenient to define Y n = Y(t n ), X n = lim t ^ 4 + x(t), Y n — lim t ^ t + y(t), 
and similarly for all higher derivatives. Of course it is X n — x(t n ) — for all n. 
With this notation, the collision rule |(6|) becomes 

(9) X n = - lim x(t). 



Next we define the variables 
(10) 



(12) 



(13) 



V>=(i/ + s-l)/2 
£ = (y-x-l)/2 

The collision condition x(t n ) = in the new variables reads 

(11) ^(**)=f(*n). 
The collision rule ((6]) becomes 

In terms of these variables the equations of motion read 

ip = -7 

i=-s-2i4 ■ 

The total mechanical energy of the idealized ball is 

£ 2 £ 2 & 

(14) E = ^ + S- + y Y +^. 
which is dissipated at the rate 

(15) f = -2^. 



The collision rule lfT2j) implies that E is a continuous function of time, even at the 
impact times. The energy is then a non-increasing function of time. 
The mechanical system has only one static equilibrium, which is 

(16) x = x = y = 0, y = 1 — 27. 
This state has energy 

(17) E m = - 7 2 /2 

which is the minimal energy of the system. 

To insure that the upper point mass at equilibrium is above the floor, we must 
require that 

(18) 7<~- 

The validity of the model may be questioned if the spring's length shrinks to 
zero, that is if, for some t, we find x(t) = y(t), or, equivalently, £(t) = —1/2. To 



avoid this event we may choose initial conditions having mechanical energy less 
than the minimum energy stored in a ball with zero length spring, that is 

1-47 

(19) E< ^^- 

Imposing an upper bound to the mechanical energy translates the fact that a real 
ball thrown on a rigid floor with excessive energy would break or undergo plastic 
deformations, thus changing the physics of the problem. However, from a mathe- 
matical standpoint, we may wish to study the case where the constraint represented 
by the floor applies only to the lower point mass. That is, we require x(t) > for 
any t, but we allow y(t) to be negative. In this paper we need not to enforce the 
restriction l|18p and (fl9| . and the only constraint is x(t) > 0. 

The main result is the following theorem, which rules out the possibility of 
inelastic collapse in our model. 

Theorem 2.1. Starting from any initial condition, there are two possible outcomes 
as t — > oo; either the lower point mass remains in contact with the floor, while 
the upper one undergoes damped harmonic oscillations; or the lower point mass 
experiences an unlimited number of instantaneous impacts with the floor, in which 
case the times of flight {5p follow the asymptotic relation 
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Although in both cases our model is well defined for any positive time, the first 
outcome is non-generic. In fact, we shall prove that 

Theorem 2.2. The contacts of the lower point mass with the floor are always 
instantaneous, that is x(t n ) = implies lim t ^ t - x(t) ^ 0, except for the solutions 
generated by a set of initial conditions having zero Lebesgue measure. Furthermore, 
this set is nowhere dense in the set of all possible initial conditions. 

A by-product of the main theorem is the following 

Theorem 2.3. Starting from any initial condition, the mechanical system described 
in this section tends to the state of static equilibrium $16\) as t — > oo. 

3. Anomalous Contacts 

The design goal of the mechanical system described in section |(2]) was modelling 
a prolonged contact between a ball and the floor with a sequence of instantaneous 
impacts, where the lower point mass reaches the floor with a non-zero velocity. 
However, there may be some anomalous contacts where this is not the case: their 
properties need to be fully understood before proving the theorems stated in the 
previous section, even if we will show that they arise only from a zero measure set 
of initial conditions. 

The first class of such contacts, that are often named grazings in the literature, 
are instantaneous contacts where there is no exchange of momentum between the 
lower point mass and the floor. A grazing event occurs when 

(20) x(t„) = x(t s ) = 0, i{t s ) > 

at some time t s . Applying the impact rule ((6]) we realize that the presence of the 
floor is irrelevant during a grazing, because in an interval around t s the dynamics 
would be the same with of without the floor. We note that a trajectory can not have 
x(t s ) = x(t s ) = 0, x(t s ) < 0, because that would violate the constraint x(t) > for 
times close to t s . 

On the other hand, it is possible to have trajectories touching the floor with zero 
velocity and acceleration, and this leads to the second class of anomalous contacts, 
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where the lower point mass maintains contact with the floor for a non-zero interval 
of time. We will call this one a sticky event. The duration of a sticky event can 
either be infinite, leading to the first of the two outcomes mentioned in Theorem 
(|2.ip . or it may be finite, then the mechanical system resumes its ordinary dynamics 
with instantaneous impacts. 

A sticky event begins at a certain time t a if 

(21) x(t a ) = x(t a ) = x(t s ) = 0. 

In such a case x(t s ) < 0. In fact, x(t s ) > is inconsistent with the constraint 
x(t) > for t approaching t s from below. Moreover, by differentiating twice the 
first of ((£]) and using the second to eliminate y and its derivatives, one finds that 
x(t s ) = implies x^ 4 '(t s ) = —7, which also is inconsistent with the constraint. 

In this situation the collision rule iJH) does not change the state of the system 
but, to avoid breaking the constraint x(t) > 0, it is clear that a new element must 
come into play, namely a non-instantaneous force exerted by the floor onto the 
lower point mass. 

Using the equations of motion |(8|) and their first derivatives, we find that condi- 
tions i|2"T|l are equivalent to 

(22) X s = X s = 0, Y s = 1 + 2 7 - 2^Y S 

where, for brevity, we define X s = x(t s ), Y s = y(t s ), etc. Moreover, x'(t s ) < 
becomes Y s < 47/1. 

The sum of the (non-dimensional) forces exerted by the spring and by gravity 
on the lower point mass is 

(23) F=±y + ny-<Y-~. 

We observe that F(t s ) = and F(t s ) = x(t s ) < 0. During a sticky event F is non- 
positive and it is balanced by the reaction exerted by the floor. The lower point 
mass maintains contact with the floor until the time t& > t s where the following 
conditions are satisfied 



(24) 



F(t d ) = 
F(t d ) > 



The second condition has to be a strict inequality because F(t) = F(t) = easily 
implies F(t) = —7, using 11251) below, and this means that the force remains non- 
positive around such a t. In the interval \t s , td] the motion of the upper point mass 
obeys the equation 

11 

(25) y + m + = -7+ 2' 

The following lemma is pivotal in order to understand the properties of anoma- 
lous contacts. 

Lemma 3.1. Sticky and grazing events are impossible if the energy of the mechan- 
ical system is less than 

3 - 2fi 2 2 

(26) E min = 2(1 + V) 7 ■ 

Proof. Using (0, J20]) or {22]), and defining z = Y s + 27 - 1, we find that the 
mechanical energy at the onset of a sticky event or a grazing is 



(27) 



Y 2 z 2 -y 2 
V ' 4 8 2 



We further have the inequality 



(28) z > 4 7 - 2pY a 

where the equality holds for sticky events and the strict inequality for grazings. 
The energy (|27l) subject to the constraint f28|) reaches its minimum for 



(29) Y s = 



4 7 M 



and 
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1 + 2/i 2 

where the energy assumes the value f26|) . For lower energies the conditions |20|) 
and l(22|) can not be satisfied, and this forbids the occurrence of an anomalous 
contact. □ 

We observe that E min > E m , the absolute minimum energy of the ball as defined 
in lfT7|) . and the coincide only in the limit fi — * 00. For all finite values of 11 there is a 
finite interval of energies corresponding to states that can not give rise to anomalous 
contacts. 

The next lemma bounds from below the duration of a sticky event. 

Lemma 3.2. The duration td — t s of a sticky event can not be made arbitrarily 
short. 

Proof. Observe that F, defined by ([23]) in terms of the solution of (|25jl. is a C°° 
function in the interval [t s ,trf]. Let us say that, with a suitable choice of Y s , the 
duration ( s — td — t s may be made arbitrarily small. Then ( s satisfies the equation 

= F(t d ) = F(t.) + F(t s )C + ±F(t s )<; 2 s + O(Cf). 

Since F(t s ) = and ( s > 0, using lj25j) and its derivatives together with (J22j) in 
order to express F(t s ) and F(t s ) in terms of Y S) we get 

2(Y S - 4 7M ) = (2 7 + Y sl i - 4 7M 2 )C S + 0(C 2 )- 

Setting a s = — Y s > we write the above expression in the form 

(30) 2a s + (2 7 - M a s )C s = 0(C). 

Then a s — > if ( s — > and both terms in the left hand side of 1(30]) are positive, 
which leads to the contradiction ( s — 0(( s ) 2 . □ 

Let us stress that an upper bound on the duration of a sticky event does not exist. 
A heuristic argument to support this statement is the following: if we take an initial 
condition where only y and y are just slightly removed from the static equilibrium 
(|2.3p . we expect the force F to be always negative, which means that the system 
evolves forever according to (f25|) . On the other hand, if we evolved backward in 
time this initial condition, an increasing amount of energy would accumulate in the 
spring, until the ball had to bounce above the floor. This means that there is a 
family of states, not in contact with the floor, that give rise to a sticky events that 
last forever. With tedious, but easy calculations, it is possible to verify that the 
state (x, x, y, y) = (0, 0, 1 + 27 — 2/ijf, y) with y — 4/x 7 — e, is the onset of such an 
infinitely long sticky event, for sufficiently small positive e. 

We can now prove the following 

Theorem 3.3. Any initial condition generates a solution which contains at most 
a finite number of sticky events. 
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Proof. Assume, by contradiction, that infinitely many sticky events occur at time 
intervals [s n ,t n ]. Clearly t n < s n+ i and Lemma (|3.2p guarantees that there is a 
positive i smaller than any t n — s n . 

The energy dissipated during the rt t,l -sticky event is 

-A* J y 2 i s ) ds - 

If |5s n | > <5 > for infinitely many n, using the boundedness of y, one sees that 
there exists < t < t , independent of n, such that \y(s)\ > 8/2 in the time interval 
[s n ,s n +t\. This implies that the system dissipates infinitely many times a fixed 
amount of energy and contradicts Lemma l|3.ip . Then Y Sn — > 0. 
We expand y 2 near s n as 

f(s n +t) = Y S 2 J + Y S Y S J 2 + ~ (fl + Y Sn Y an ) t 3 + (t 4 ) 

and we observe that lf25|) and l(22|) yield Y Sn — ~2~f. Setting e„ = Y Sn it follows 
that there exists < t < i and independent of n such that 

y 2 (s„+t) > at 3 - Ce n t 

in [s n ,s n + i\, for suitable a, C > (any a < I6/37 2 works). Then the system 
dissipates at least bi 4 — De n P, b — a/4, D = C/2, during the n th sticky event. 
Since e„ — » 0, if n is big enough this quantity is bounded from below and, as 
before, we are in contradiction with Lemma (|3.ip . □ 

Theorem 3.4. Any initial condition generates a solution which contains at most 
a finite number of grazings. 

The proof of this theorem is trivial, once the truth of l|2.ip and l|2.3p has been 
proven. For coherence with the topic of this section, we give it here, rather than 
at the end of section \E[ Of course, this theorem is not used in the following two 
sections. 

Proof. The absence of inelastic collapse (section HJ forbids an infinite sequence of 
grazings taking place in a finite time. By Theorem (|2.3p the asymptotic state is the 
static equilibrium lfT6|) to which corresponds an energy smaller than the minimal 
energy (J26j) needed for a graze. Then an infinite sequence of grazings in an infinite 
time is also forbidden. □ 



4. Absence of Inelastic Collapse 

4.1. Preliminary Considerations. In this section we prove that the idealized 
ball described in section [2] does not experience infinite contacts with the floor in a 
finite time. The proof focuses on true impacts, that is contacts such that X n ^ 0, 
so we must briefly discuss the case X n = 0. 

As we have already seen, this leads to anomalous contacts, either grazings or 
sticky events. In the latter case, if the duration of the event is infinite, the problem 
of inelastic collapse is avoided, because equation (|25| holds forever. Otherwise, the 
discussion of this section may be thought as taking place after the last sticky event, 
thanks to Theorem (|3.3p which allows only a finite number of such events. 

Grazings, too, may not happen after an arbitrarily high number of contacts, but 
at this stage there isn't a simple and quick way to prove this statement. Instead, 
for each of the cases considered below we will take care to show that X n > for n 
large enough. 
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4.2. Asymptotic relationships at collapse. In order to show that there is no 
inelastic collapse, we assume, by contradiction, that t^ is a finite positive number. 
This assumption leads to explicit approximate expressions for the function y and its 
derivative for t close to too- In turn, these expressions allow to deduce asymptotic 
laws for X n , X n , X n , and r„ which lead to contradictions. 

From (jT5j) we have that the mechanical energy of the system is a non-increasing 
function of time, and, taking into account that it can never decrease below the 
energy l[T7j) of the static equilibrium (fTo) we conclude that the mechanical energy 
has a finite limit. The existence of such a limit implies that the dynamical variables 
tp, tp, £ , £, appearing in the expression i|14p for E, are bounded, and so are their 
linear combinations x, y, x, y. 

The equations of motion ([8]) guarantee that x and y are also bounded. In general, 
iterated differentiations of (HJ) show there exist positive constants Ck such that 



(31) xW(t) , yW(t) 



,0)i 



for every t € U n (t n ,t n+1 ). 

By construction, the positions x and y of the point masses are continuous func- 
tions of time. The collision rule (J6]) makes x discontinuous at the impact times t n , 
while the rule (0 states that y is continuous. Then x £ C° [0, too) and y € C 1 [0, too). 
Furthermore, the limits Yoo = lim 4 _> too y(t) and Yoc — licit— ^ y(t) exist and are 
finite, because y <G C 1 [0,t oo ) and y is bounded. 

Now we prove that position and velocity of the lower point mass go to zero as 
t — > too ■ Integrating the first equation in (fT3|) and using (HH) we find 

(32) 7 r„ = - E n+1 
which may be rewritten as 

(33) 2 7 r„ + Y n+1 - Y n = X n + X n+1 . 

Because the l.h.s. tends to zero as n — > oo, and X n > 0, then we have limn^oo X n — 
0. But, considering that X^ are bounded quantities, for t € [t n , t n+ i) we may write 

x(t) = X n (t - t n ) + 0{tI) 

and 

x(t) = X n + 0(r„). 

Therefore we have: 

(34) t lim ±{t) = lim x(t) = 0. 

Furthermore, recalling that y is bounded, we have that Y n+ i — Y n — 0(r n ). Ob- 
serving that X„ is never negative, from 1)331) follows 



(35) X„ = 0(r n ). 

Having assumed that J2 n r « converges, we must conclude that J2n X n converges as 
well. 

(k) (k) 

The expressions for Xn and Y„ , found by differentiating repeatedly the equa- 
tions of motion jSJ, may all be written as continuous functions of X n , Y n , Y n , which 
have a limit. Then there also exist lim„_ ) . o0 X n k ^ = X^ and lim„^oo Y^ — Y^\ 

Higher derivatives of x and y are generally discontinuous in t„, but their jumps 
are all amenable to the jumps of x. To prove this fact we observe that the collision 
rule J6]) may be written as 

(36) [£(*„)] = 2X n 
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where [£(£„)] = x(t+) — x(t~) is the jump of the function x around t„. Then, from 
the first of equations J8|), recalling that x € C°[0, too) and y € C 1 [0, too), it follows 
that 

[x(t n )} = -2fiX n 

and from x + y = —2~f we have 

[y(t„)] = 2^X n . 

Analogously, differentiating equations with respect to time, and using x'+y = 
we may evaluate the jumps of the third derivatives: 

[x(t n )]=-[y{t n )] = (4 t i 2 -l)X n . 

Iterating the procedure we find that there exists a sequence of real numbers c k such 
that 



(37) 



? (fe) (*n)] = -[y (fc) (i™)] = c k x 71 



Next we observe that 

oo 



h—n 



,(k+ 



h=n+l 



Using (J3TJ and J37J) we have 



(38) 



oo n 



< c k+1 T n + c k V n 



where we have defined the time left to collapse 
(39) T n — t n 

and the residual sum 



OC 

V fe+1) (*)^ + J2 [y (fe) (^)] 



(40) y„= ^ x h . 

h=n+l 

Taking t n <t < t n+ i, from l|3Tjl and ([381 evaluated at n + 1, it follows 
(41) 

Observing that 
and using (|4lj) for fc + 1 we have 

*1* } -V (fc) (*) = (>l fc+1) + O(ioo - a)) <fa+0(V n ) = Y^+V^-V+O ((*„, - t) 2 )+0(K) 
hence 

(42) yW (t) = y« + ll fc+1 )(i - too) + 0((t - too) 2 ) + O(K). 

Choosing fc = 2, by repeated integration of l(4*2|) . exploiting the continuity of y and 
y up to too, we obtain 

(43) y{t) = too + Yocit - toe) + 2 Yoo(i - ioo) 2 + 0((t - too) 3 ) + 0((*oo - t)V n ) 
and 

(44) y(t) = Foo + too(t - too) + 2^oo(< - ^oo) 2 + 

+ iFoo(t - too) 3 + 0((t - too) 4 ) + O((too - tfV n ) 


10 



Using lf43|) and I144II in the first of equations (J8j we find that the motion of the 
lower point mass, under the assumption of inelastic collapse, would be described 
by the following equation 

(45) i + fix + -x = P - -t) + R{too - tf + 0((too - t) 3 ) + 0((too - t)V„) 
where 

P 

(46) Q 

R 

Because P, Q, R are numbers, the motion of the lower point mass is now decoupled 
from that of the upper point mass, and equation (|45| is closed. Furthermore, we 
may evaluate the equations of motion ((8J and their first two time derivatives as 
i — ► i+. Using (03l), (05) and (02) for k = 2,3 in order to eliminate y and its 
derivatives, observing that X n = 0(T n ) because of <f35|l . we find the following 
asymptotic relationships 

( X n = P- ixX n - QT n + RTl + 0(T 3 ) + 0(T n V n ) 

(47) I X n = Q - fiP + (m<3 - 2R)T n + 0{Tl) + 0(V n ) . 
[ X^ = (fi 2 - \)P -fiQ + 2R + 0(T n ) + 0(V n ) 

4.3. Absence of Inelastic Collapse. In the interval between two consecutive 
impacts, the motion of the point masses is described by smooth functions. Thereby, 
there exist two numbers Cn and r] n , with t n < £ n , rj n < t n+ i, such that the following 
expressions hold for t n < t < t n+ i 

2 3 

(48) x{t n+l ) - X n r n + X n ^ + 
and 
(49) 

Since x(t n+1 ) 

(50) X n r n = -2X n - ic-(Cn)y- 

Inserting this expression in (09]) and applying the collision rule ([9]), we obtain the 
following map for the speed of the rebounds 

(51) X n+1 =X n - A„t-2 
where 

(52) A n = {x{r) n )/2 - x(Cn)/3) . 

With the help of the map (f5lj) . we shall examine the dynamics generated by 
equation (05| separately for the three cases P ^ 0; P — 0, Q ^ 0; and P = 0, 
Q = 0; that is i„ ^ 0; = 0, Xoo ^ 0; and i"^ = 0, Xoo = 0. They 
are exhaustive of all possible asymptotic behaviors of the ball because, from the 
definitions (06j, P = 0, Q = implies R = —7/2. We shall see that they all are in 
contradiction with the hypothesis of inelastic collapse. 
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y 



T 

x(tn+i) = X n + X n r n + x(r] n ) — , 
= 0, from (081) we find 



4.3.1. The case Xoo ^ 0. Recalling that x is a bounded function, and that lim, woo t„ 
0, dUJ yields 

(53) r„ = 0(X„) 

and also 

X 00 = lim [-2— + 0(t*)| <0. 

n—oo y Tn J 

This implies X n > 0, X n < and r„ ~ X n for large n. Then (|51"Tl gives 

X n+1 = X n + 0{Xl). 

From Lemma 1|9.3|1 follows that ^2 n X n = +oo, then J2 n Tn diverges as well. Thus 
the inelastic collapse can not occur in this case. 

4.3.2. The case X^ = and Xoo ^ 0. Equations (@8} through (|52f. are still 
valid because they do not depend on the asymptotic behavior of X n and X„ ■ In 
particular, in the limit n — > oo, since 'i'(Cn)> x{r]n) = Xn + 0(r n ) in the definition 
(|52"ll . we have A„ — > Xoo/6. If we assume Xoo < 0, the map l[51"j) does not allow 
X„ — > 0, which is inconsistent with (f50|) and r„ — > 0. Then Xoo 7^ implies jf„ > 
for large n. Furthermore, A„ > 0, together with r„ > implies X n > 0. In fact, 
if X n = for some n, the map (f5lj) would give X„+i < which is in contradiction 
with the constraint x(t) > imposed by the presence of the floor. 

Using X n ,T n ,X n > 0, from equation 1(50]). we find that 



In < 

for all n large enough. Furthermore, the hypothesis X^ = implies 

(54) lim — = 0. 
Observing that T„ = T « we a l so have 

(55) hm ^ = 0. 

n— >oo ± n 

Using l|5ip and (|39| we can write the following map for the pair (X n ,T n ): 

X n +i = X n — A„t^ 

In order to have an inelastic collapse we need to find at least one sequence of pairs 
(X n , T n ) — > (0, 0) as n — > which satisfies this this map and l(55|). 

Having proved that X n ,X n > and X n < for n large enough, recalling that 
x{t n+ \) = 0, the time of flight computed from ()48| is 



(56) 



(57) T _ 3 \ X n\ /\ L 8x(( n )X n 



2 x (C„) ^ V 3 X 
With ([56]) and (J5TJ) we compute 



(58) ^~5" +1 =4- A ' 



x„ T n — r n+ i 



8 ^nKnlXg 

3 XI 



From (|4*7|l and (|55|) . recalling that P = Xoo = 0, we have X„ ~ —QT n and X„ ~ Q, 
that is A„ ~ Q/6. Using these asymptotic expressions in (f58|) we find 

/ \ T n X n — X^j-i 2 
(59) hmsup-^- <- 
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which, in turn, implies 



(60) £±i > £ 



This is in contradiction with (|55|l . 

To conclude the study of this case we observe that the second non-zero root of 
(|4"5f. namely 




(61) fB = rj^ii + ./i_ 8 ' if ^)- Y « 



2 



2 af(Cn) V V 3 ^ 



can not be the time of flight for large n, because X„ <~ —QT n and X n ~ Q would 
lead to the absurd f n > T„. 



4.3.3. T/ie case X^ — and Xoo — 0. In this case we have P = 0, Q = in ((33 
so that jf„ = o(T n ). Looking at (|3H)l now we find 

(62) X„ = o(r„T„) 
which implies that the residual sum $40$) is 

(63) V n = o(T„ 2 ). 



From (|46| . with the help of the equations of motion ((8J and their derivatives, we 
have R = —7/2. Using these relations back in fTfj) we find 

( X n = -i-yT^ + o(T^ 
(64) \ X n =lT n + 0{T*) . 

{ Xi 4) =- 7 + 0(T„) 

The first of these equalities show that X n < for large n. Then X n > for equally 
large n. In fact, X n — and X n < imply x(t) < for t close to t n , which is 
incompatible with the floor constraint x(t) > 0. 

Evaluating x(t n+ i) with a power series expansion around t n , we have 




(65) V „=-(^r„ ^ + ^^+0(7-*). 

In the same way, a power series for x(t n +i) together with the collision rule Q gives 

(66) X n+1 = -X n - X n r n - \x n r 2 n - \xWt* + 0(r*). 



We use l(65|) in f66|) in order to eliminate X n and X n+ \, we substitute the higher 
derivatives with the expressions ([64]) . and we write r n = a„T„ for some a„ G (0, 1). 
After some algebra, this yields the following implicit map 

(67) (1 - a n f (6a„+i - Aa 2 n+l + a£ +1 ) = 6a„ - 8a 2 + 3a s n + o(l) 

where the term o(l) goes to zero as n — > 00. We prove in Lemma (|9.4p that if the 
sequence {«„} stays in (0,1), then it converges to zero, that is r„ = o(T n ). Of 
course, a n > 1 for some n would lead to r„ > T n , which is an absurd. 

Having proved that, for large n, it is X n > 0, X n < 0, X n > 0, it follows that 
expressions l|55p through (|58|) still hold in this case. Using (|64|) in ([65I we find 
~ jT 2 r n . Then, recalling that 'i'(Cn) = X„ + 0(r n ), we find 

^ (Cn)^ n T n 
13 



and we conclude that the square root in ([58| approaches one for large n. Using 
A„ ~ J„/6 and ([6^1 in ([58| . in this case gives 



lim 

n — >oo 



Thus we have again 



T n X n — X n+ i _ 2 
X n T n — T n+ \ 3 

X n +i X n 
T n +i T n 



which contradicts (|33)l . Finally, we observe that, also in this case, ([6"T]) gives r„ > 

5. Asymptotic State 

We now carry on to prove Theorem l|2.3p . If the asymptotic dynamics is sticky, 
then equation <f25|l holds for all times greater than the last attachment time t s . It is 
straightforward to verify that as t — > +oo the system reaches the equilibrium (fl6|) . 

If the asymptotic dynamics is not sticky, then a rigorous proof is not trivial. A 
useful information, more or less implicit in the calculations of the previous section, 
is that 



(68) 



X n = 0(T n ^) 



To see that this equality holds, we remark that there exists a number rfo-i € 
(i„_i,t„) such that 

-X n = A„_i + x(?7„_i)r„_i 

which implies l(68|) . since X n is not negative for any n. 

We split the proof in two lemmas, and a final theorem. Both lemmas, ideally, 
are nothing more than a Lyapunov fixed point theorem, but they have to deal with 
technicalities arising from the presence of the impacts. 

Lemma 5.1. // the asymptotic dynamics is not sticky, then lim^oo £(t) = 0. 
Proof. We assume by contradiction that 



(69) 



lim sup 

t — >oo 



> 0. 



If ([69]) holds, having excluded the presence of inelastic collapse, it is possible to 
find a monotonically growing sequence of times Sk — * oo, with Sk + i — Sk > 2 and 



limsupfe^ £(s fe ; 

recalling that E(t) 
such that 

(70) 

and 

(71) 



> 0. Because £ and all its derivatives are bounded functions, 
= -2ni 2 (see eq. HDJ we see that there exist positive S and C 

E(s k ) < -5 



E(t) 



< C 



for all k and for all t S U^L (^n> t n +i)- We claim that there is a t such that 
|* — Sfc > 5/(16C) and E(t) < —8/2 for any t that is not a time of impact in 
the interval having t and Sk as extremes. Let t„, t n+1 be the times of consecutive 
impacts that bracket su'- to prove the claim we have to consider two different cases: 
1. If t n+ i — t n > 5/ (8C), then there exists t S [t n , such that 

(72) \t-s k \ = 
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16C 



From (f7l|) . for any t between t and Sk we have 

E{t)-E{s k)=H 

t - Sk 

for some £ between t and Sk- Observing that \t — Sk\ < <5/(16C), from (|70|) we get 

E(t) < E(s k ) + — < -- 

2. If 

(73) t„+i - i n < — 

then let m be the integer such that 

( 7 ,\ \ tn+m+l — t n > 5/ (4C) 

K ' 1 t n+m -t n < 6/ (AC) ■ 

There exists t £ [t n + mi t n+rn+ i] such that 

(75) f -*" = ^' 
Recalling that t n < s k < t n+1 then 

t — Sk > t — t n+ i = (t — t n ) — (t n+ i — t n ) 
and from ((73]) and ((74j) we get 

(76) f " Sfe> ^7 

From lfT0|) we have £ 2 = (i 2 — 2xy + y 2 ) /4 and from ([36| follows that the jump of 
£ 2 across an impact is 

(77) [i 2 (t n )} = -X n y(t n ). 

Recalling that E(t) = — 2^ 2 and that X n = 0(r n -i), we have can choose the 
constant C in (|7T|) in such a way that 

(78) [E(t n )} < Cr n _ a 

for all n. Moreover, for t € [s k ,t\ such that t n + s < t < t n+s+ i for some s <m, 

(79) E(t) - E(s k ) = (E(t) - E(t+ +S j) + (E(t- +S ) - E(t+ +S _ 1 )) +■■■ 

S 

■■■+ (E(t- +1 ) - E(s k )) + Y^[E(t n+i )}. 

Using (J7TJ) and (JTHJ) we obtain 



i=i 



£(t) - S(s fe ) < C (t - s k ) + C (t n+m - t n ) < 2C (t - t n ) = - 



and, from ([TOJ) 

E(t)<- 6 -. 

The claim is proved, and for any k there exists an interval not shorter than 
<S/(16C) having Sk as an extremum. During each of those intervals the power 
dissipated is at least —8/2, which gives an infinite dissipation of the mechanical 
energy, in contradiction with the fact that E(t) has a lower bound. Then ([69)) is 
false, and we must conclude that lim^oo £(t) =0. □ 

Lemma 5.2. // the asymptotic dynamics is not sticky, then lim^oo^) = —7. 



Proof. The proof is identical to that of the previous lemma, except that, in place 
of the mechanical energy, we use an ad hoc Lyapunov function L defined as 

m = \[(£+ + i)f + $ + + t)) 2 + e\ + 

Because ip, £ and all their derivatives are bounded, £ and all its derivatives are 
bounded as well. We observe that £ is a continuous function of time, even across 
the times of impact. Furthermore we find 

(80) £(t) = -2M£W+7) 2 + O(0. 

We assume by contradiction that limsup^^ +7I > 0. Proceeding as in the 
previous lemma, and exploiting the fact that lim^o £(t) = makes f80|) negative if 
|£(t) — 7I is bounded from below from zero, we conclude that lim^oo C(t) = —00, 
which contradicts the boundedness of C. Therefore lim^oo (£(f) + 7) = 0. □ 

5.1. Proof of Theorem (12.31) . From the collision condition (JTTJ) and lemma (|5.2p 
it follows 

lim ip{t n ) = -7. 

n — >oo 

From the collision rule (TT2l) and lemma (15.11) it follows 



lim = 

n — >oc 

and from ([H]) we deduce that limn^oo E(t+) — —7/2. Since the energy is a contin- 
uous, decreasing function, then 

lim E(t) = -1 

t — >oo Z 

The only state with this energy is the state of static equilibrium (fT6| . 

5.2. Proof of the Main Theorem pTT) . 

Proof. Theorem (f273|l . together with {32]), yields 
(81) lim r n = 

n — '00 

By Theorem l|2.3p we have also 

lim y(t) = 1 — 27, lim x(t) = lim x{t) = lim y(t) = 

t — >00 t — >00 t — >OG t — >00 

From the equations of motion (6) and their derivatives we can easily compute 
lim X n = -27 , lim X n = 2^7 

n — ^00 n — >oo 



We follow here the steps of section 14731 in order to write a map for X n . From l(48|) 
we get t„ as in lf57j) (the other non-zero root of l[48|) would contradict ([81])). We 
note that for £„ e (t n ,t n+1 ), x (( n ) =X n +0(r n ), so 

(82) lim x (C„) = 2//7 

n — >oo 

Using l|49p and (|50]l we can write 

(83) X n+1 =X n - A n r^ 
with A„ defined in ([52]) . From (|82]l we have 

(84) lim A n = ^ 

n — >oo o 

Then, expanding the square- root in ([57]) . from (|83]l we get the map 

(85) X n+1 =X n -a n Xl + 0{Xl) 
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with a n — 4 A„/ X%- Since lim^oo a n = /i/(37) we can rewrite ((8 
(86) X n+1 =(l- f-X n ) X„ + »( X~ i 

Applying Lemma (|9,ip to <f86|l we get 



if 3 ^ 1 
/Lt n 

Dividing both sides of (f5Q]) by r n we get 



then 



X n Xn , ■■■ / j. \ 

T~n 2 



lim = 7 



which proves the thesis. □ 
6. Non-Genericity of Zero- Velocity Contacts 



In this section we give a proof of Theorem 112.21) . The state of our idealized ball 
is a vector u S M 4 . Without loss of generality, we may assume that the first two 
components of u are the position and the velocity of the lower point mass. We 
define the hyperplane C of contact configurations as 

(87) C= {ueR 4 |u-ei =0} 

and the plane S of contact configurations with zero velocity 

(88) S = {ueC|u-e 2 = 0}. 

We observe that the vector field \ that defines the equations of motion JH|) , on 
the hyperplane C is 



-[IX - 7 - ^ + 2V + M 

y 

V -m - \v - 7 + \ + j 



It follows that x(u) is transverse to C if and only if the second component of u is 
non zero. That is, \ ls transverse to C everywhere but in S. 

To prove that the set of initial conditions generating orbits that reach S in a 
finite time has zero measure, we need to define the backward return map of the 
flow on C. This requires some care because the points of S are contiguous to points 
that cannot be reached by the dynamics. 

The time evolution of state vectors in the interval (t n ,t n+ i) is determined by 
the flow <p : R x M 4 — ► R 4 of the equations of motion J8]) . We may define an inverse 
time-of-flight function implicitly as the largest negative number r~such that 

(89) i^(r",u)-ei = 

for any ueC", where C~ is the subset of vectors of C with x < 0. Some care is 
needed in order to define r~ for the states belonging to S. We observe that if x = 
and x > we have a grazing event, and definition (|89| is applicable. If 'x = then 
we need to evaluate the force (f23| exherted by the ball on the floor: if F(u) = 
and F(u) < then the state u is the beginning of a sticky event, and definition 
([89]) is again applicable; for all other states on S a sticky event has already begun, 
and we define t~ = 0. 

We observe that the backward-in-time evolution of a sticky state does not leave S 
as long as t~ = 0. On the other hand, any state with t~ ^ leaves C following the 
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u_ 2 =<p(T:_ n Ru_,) 





u_,= (p(T ,u ) 



x=0 



X 



R 



Figure 2. In proximity of a non-sticky state uo that evolves back- 
ward in time into a state u_i ^ S, a subset of S (the thick segment) 
is mapped into a subset of C of zero measure (the thick wiggly line) . 
Furthermore, if u_i evolves backward into u_ 2 S then is a 
diffeomorphism between a neighborhood of u_i and a neighbor- 
hood of U_2. 



flow f backward in time. Then, for these states we introduce the transformation 
T- :C+ -> C+ defined as 



where R is the reflection matrix that changes sign to the second component of u. 
A sketch of this process is given in figure El 

Let us assume that a given pair (t ~ , u ) satisfies l(89| with ^ 0. If v?(tq" , u ) ^ 
S, then 



and the implicit function theorem yields a smooth mapping r~ in an open neigh- 
borhood of (t ~, u ). However this map need not coincide with our inverse time of 
flight in the whole neighborhood, since it depends on the free flow tp which does 
not take into account the presence of the floor. To overcome the problem, we ar- 
gue as follows. If Uo £ S we may not apply the transformation ([90)1 to an entire 
open neighborhood of Uo in C, because that neighborhood would necessarily contain 
points of C~ , where ([90| is not defined. Then, if Uo is a grazing, we consider a small 
open neighborhood U C S of u in S containing only grazing states and satisfying 
T~ (U)C\S = 0. Since Ru — u if u g S, we deduce that the restriction of T~ on U is 
as smooth as ip. If Uo is the beginning of a sticky event, then we can not consider an 
open neighborhood of Uo in S because it would contain states having t~ = 0, and 
for those the transformation l|90p is not defined. Then we consider open intervals 
U along the straight line defined by x = 0, F = 0. If U is small enough, then it 
does not contain states for which F > and moreover T~(U) n S = 0. Then the 
restriction of T~ on U is as smooth as ip. 
The union of all the above sets U, U is 



The set S is contained in a two-dimensional plane and has zero measure in C + , which 
is three-dimensional, and the previous discussion proves that T~ is as smooth as <p 
on S. Then the image of S in C + through T~, namely T~(S), is also a set of zero 
measure in C + , see e.g. p~4l Chapter 2, Proposition 1.6]. By analogy, for n > 2, we 

18 



(90) 



T-{u)=tp(r-,Ru) 




a 



S= {ueS | T"(u) is defined and T"(u) £ S} . 



define 

5" = {ue T-iS 11 - 1 ) |T-(u) ^5} 

and we set S 1 = S, S° — S. As before, T~ is smooth on <S" and then, by induction, 
S n+1 C T~(S n ) has zero measure in C + for every n. By the sigma-additivity of 
the Lebesgue measure, it follows that 

£ = u~ 0< S" 

has zero measure, too. 

The bouncing ball starting from an initial condition u £ R 4 , with u • ei > 0, will 
experience a grazing or a sticky event if its first contact state v?(r(u),u) satisfies 

¥>(t(u),u) e E. 

So, the set A of anomalous initial conditions that reach S after a finite number of 
contacts is contained in 

U t < <^,£) = <p((-oo,0] x S) 

which has zero measure in R 4 , by the smoothness of ip with respect to all its 
variables. 

We now prove that A is nowhere dense in R 4 . Because A has zero measure, 
no sphere centered in a point belonging to A is a subset of A. Then we need to 
show that around each non-anomalous points there is a neighborhood that does not 
intersect A. 

Let us assume that u ^ A . Let us call u(i) the state of the system at the time 
t, starting from the initial condition u. From the results of section [3] it follows that 
after iV contacts it is i£(u(ijv)) < E min , where E min is the minimal energy required 
to have an anomalous contact (see Lemma (|3.1[) \ Since the mechanical energy is a 
continous function of the state of the system, there is a neighborhood V of u(t/v) 
such that V n S = and E(v) < E m i n if v € V. Furthermore, V may be taken 
small enough to ensure that T _n (V) fl S = for n = 1, • • • , N. Finally, we observe 
that any neighborhood U of u such that ip(t(U),U) C T~ n (V) does not intersect 
A. 

7. Numerical Results 

7.1. On the Nature of the Restitution Coefficient. The dynamics of the 
mechanical system described in section (J2]) may be easily simulated on a digital 
computer, finding the instants of impact of the lower point mass with the floor 
with an accuracy as high as the machine precision. 

If the mechanical system is left free to fall with its internal degree of freedom 
not excited (that is with initial conditions £(0) = £(0) = 0, "0(0) > 0, and arbitrary 
4>(0)) the resulting dynamics may, initially, approximate that of the rigid ball model 
with constant restitution coefficient (eq. {J}), provided that the spring is weakly 
dissipative and sufficiently rigid, and that the energy of the initial condition is 
sufficiently high. An example is illustrated in figure (J3]), where the mechanical 
system is left free to fall from an height of six non-dimensional units. 

This fact is explained by observing that, upon an impact, the vibrational and 
translational velocities exchange with each other (eq. l(T2)l ). so the center of mass 
of the ball has zero velocity just after the impact. In this situation we recognize 
three characteristic times 

T^, is the free-fall time of the center of mass, T£ is half the proper period of the 
spring, and Td is the e-folding time associated with energy dissipation from the 
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Time 

Figure 3. Numerical simulation with initial conditions £(0) = 
£(0) = ^(0) = 0, V(0) = 6, and 7 = fx = 0.01. The upper curve 
is y(t), the lower curve is x(t). The solid dots show the instants 
when there is an impact. The inset magnifies the region within the 
dashed rectangle. 



spring. If Tjf, is much larger than T% then the center of mass does not have enough 
time to gain any appreciable speed before the expanding spring causes a second 
impact, which exchanges vibrational and translational velocities again. After this 
second impact the mechanical system is again propelled upward, with a very small 
vibrational energy left in the spring. At this point, if the kinetic energy of the 
center of mass is high enough to cause a time of flight to be of the order of Td or 
larger, then any residual vibrational energy will be dissipated to negligible levels 
during the flight, and the mechanical system will present itself to a third impact 
with the spring essential at rest, causing the double-impact dynamics described 
here to repeat itself. 

In the limit 7^0 the residual vibrational energy tends to zero, and all dissi- 
pation happens during the double impact. Then it is appropriate to introduce the 
restitution coefficient 

(92) r = *„/*„ +2 

where ^f n and ^l n +2 are, respectively, the velocity of the center of mass before and 
after one of these double-impact events. It can be shown (see [15] for details) that 
it is 

(93) r = e-^. 

If we take our mechanical system as model of a bouncing ball, we are lead to 
say that during a double-impact event the ball is in contact with the floor: it takes 
this finite amount of time for the center of mass of the ball to exchange momentum 
with the floor and revert upward its velocity. Only flights longer than are to be 
taken as macroscopic flights. In the regime discussed here, our model may be seen 
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as a version of the rigid ball model of eq. {TJ in which the contacts have a finite, 
rather than infinitesimal, duration. 

When, after a number of double- impact events, the time of flight becomes shorter 
than the characteristic dissipation time Tj, subsequent impacts will be able to store 
a significant amount of energy into the vibrational mode, subtracting it from the 
energy of the center of mass. This can be seen clearly in figure ([3]) where there is 
a sudden drop of the maximum height reached by the mechanical system between 
the fourth and the sixth macroscopic flight, and, at the same time, the trajectories 
of the two point masses become rather complicated and far from parabolic. 

It is interesting a comparison between experiments in which a spherical bead 
is free to bounce repeatedly upon a rigid surface until it comes to rest, and our 
numerical simulations. Because the velocities of the center of mass of a bead are 
difficult to measure, it is customary in this kind of experiment to measure the time 
of flight, which is more accessible to the observer. In the experiments reported in 
|13j the coefficient of restitution is expressed as 

/ n .\ fn+l 

(94) r n = — — 

Jn 

where /„ is the time of flight of the bead, (which does not include the time of 
contact of the bead with the rigid surface). If the rigid ball model of eq. ([I]) is 
valid, then f94|) is equivalent to the standard definition of restitution coefficient as 
ratio of velocities before and after the impact. If the bead is not rigid, then l(9i|) 
must be taken as an alternate definition of restitution coefficient. In our numerical 
simulations, taking into account the interpretation of double-impacts as contact 
time with the floor, it is natural to define 

(95) /„ = T 2n . 

To simulate the very rigid tungsten carbide bead used in [13] we set 7 = 10~ 5 , while 
(j, = 0.01 leads, according to the approximation lj93j) . to the restitution coefficient 
r = 0.969 • • • , which is about the same as what is measured in [13] for relatively 
high- velocity impacts. The simulated mechanical system, as the experimental bead, 
is left free-falling from a height equal to one quarter of its length. The restitution 
coefficient computed from the simulation, according to the definition ([9^1 . is shown 
in figure (01). Initially the value of r is constant, and it is in excellent agreement with 
(|93| . Then r declines, and finally it begins to fluctuate without a clear pattern. At 
even later times (not pictured in figure) the double-impact dynamics is completely 
disrupted, as in the case of figure J3}, and the definition (f95|) ceases to be meaningful. 
This kind of behavior is qualitatively the same as that of the bead observed in [13] 
(see their figure (7)). In their case, after fluctuations in the value of r of the 
same order of magnitude as those of figure (@]), the bead is observed to vibrate 
while maintaining contact with the floor. We stress that the behavior of figure 
(|4]) remains qualitatively the same for any choice of parameters corresponding to 
slightly damped, very rigid springs. 

A different way to look at the same dynamics is pictured in figure where we 
show the sequence of states belonging to the set C of contact configurations |87|) 
for the same simulation that generated figure Q . Initially, the dynamics is clearly 
dominated by double-impacts, whose signature is the alternation of negative and 
positive values of y, which gives a characteristic zig-zag look at the early part of the 
sequence. At later times the sequence of states becomes very disordered, and lacks 
any easily recognizable pattern, except a tendency to move towards the equilibrium 
point. Figure ([6]) is a magnification of figure |(5j) close to the equilibrium point. 
At these even later stages, the dynamics is well approximated by the map (|85| . 
The length of the spring and the velocity of the upper point mass are both subject 
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Figure 4. Restitution coefficient computed using the definition 
f94|) from a simulation with 7 = 10~ 5 , /1 = 1CP 2 , and initial con- 
ditions £(0) = £(0) = ^(0) = 0, V(0) = 0.25. 




Figure 5. The dots represent the state of the system immediately 
after a contact with the floor (that is, x = 0) for the simulation 
of figure ((4J. The equilibrium position of the upper point mass 
is y eq = 1 — 27. The thin lines have the sole purpose of making 
evident the temporal sequence of states. 
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Figure 6. Magnification of the figure (HJ, close to (0,0,0). 




Figure 7. Same as figure ([5]), but with (i = 2. 



to damped oscillations around the equilibrium position, and this gives a spiraling 
appearance to the sequence of states. 

The reader may have noticed that the asymptotics of Theorem l|2.ip are indepen- 
dent of the particular value of the damping parameter /i. In particular, the results 
of the theorem apply equally to under-damped and over-damped oscillations. In 
figures J7]) and (JSj) we show a simulation that uses 7 = 10~ 5 , /1 = 2 and the same ini- 
tial conditions as the simulation of figure (J5j) . The length of the overdamped spring 

23 




10 5 ± 



Figure 8. Magnification of figure 0, close to (0,0,0). 



shrinks for the first few impacts, then slowly expands, approaching the equilibrium 
value monotonically from below. The system enters immediately the asymptotic 
regime, and the is no disordered transient. From a macroscopic point of view, the 
idealized ball may be taken as performing a single, totally inelastic impact with the 
floor. 

7.2. Sticky Solutions. We have performed several numerical simulations using 
the following initial conditions 



For e = they obey the conditions l(22|) . leading to a sticky event which lasts for 
t c = 5.23 • • • non-dimensional time units. In figure |(9|) some sample trajectories 
are shown. It appears that, as e decreases, the solutions tend to the solution with 
e = 0. To quantify the convergence, we have computed the norm 



where y e and y s are, respectively, the position of the upper point mass with e ^ 0, 
and with e = 0. The results are shown in figure (JTOj) - The root-mean-squared 
difference IpT)) falls to zero faster than linearly with e. 

The numerical evidence supports the idea that the dynamics during a sticky event 
may be approximated up to any degree of accuracy by a non-sticky dynamics. 



Our simple model shows that taking into account explicitly the deformability 
of a bouncing body avoids the pathology of the inelastic collapse. Of course this 
comes at a price: our description of the bouncing ball, in spite of its simplicity, is 



(96) 




1 + 2 7 - 2^Y 
-0.1 



(97) 




8. Discussion and Conclusions 
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Figure 9. Simulations using the initial conditions ([96)1 for three 
different values of e. Solid dots ifiark the times of impact. In the 
panel e = 0, in the interval between the first two solid dots, the 
lower point mass is in contact with the floor. 
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Figure 10. Number of impacts performed by the lower point mass 
in the interval (0, t c ) and root mean squared difference (|97|) as 
function of e. 



considerably more complicated than a restitution-coefficient based model. This is 
affordable for one, or maybe a few, particles, but it is unbearable for systems with 
a very large number of bodies, even resorting to numerical simulation. 

One may be tempted to overcome this difficulty by returning to a description 
based upon a restitution coefficient of the form 

(98) r(u) = [ 6XP ' U>UC 

{ 1 - £jU, u<u c 

where u c = 37/i~ 1 ^1 — exp(— pmj^J\ — ^ 2 )^j . This is obtained by stitching the 

expression |93|) . which is valid for times of flight larger than the characteristic 
damping time of internal vibrations, and the coefficient in expression l(86|) . which is 
the map that approximates the asymptotic vibrations of a ball close to equilibrium. 

Although (|98| is very accurate for large and for small values of u, a glance at the 
figures of section [7] shows that (for underdamped springs) there is an intermediate 
range of impact speeds where f98|) simply fails. When the internal vibrational 
mode is excited, but the times of flight are not yet infinitesimal, the dynamics is 
too irregular to be parameterized with a velocity-dependent restitution coefficient. 

We may hint that a satisfactory parameterization will need some form of book- 
keeping of the energy stored in the internal vibrations, and some rule to determine 
the exchanges between that and the energy of the center of mass. At this stage, 
we are unable to suggest anything more specific than these vague considerations. 
In particular, it is not at all clear if accounting for internal vibrations will help 
the ongoing search for some sort of hydrodynamic limit for systems with a large 
number of dissipative particles. 

We remark that the behavior of the bouncing ball, close to equilibrium, does not 
depend crucially on the linearity of the spring. A nonlinear variant of the model 
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described in section <[2j) is embodied by the equations 



i> = -7 

e = - P ^r-2^i£i fc 

with p, n > 0, a, 6 > 0. The impact times t n are determined by condition (fTTjl . 
where the collision rule 11121) must be applied. The constraint tp(t) > £(t) holds 
at all times. We recover the linear model for a = b = 0. It has been argued, on 
the basis of Hertz's contact law, that the correct choice for modeling homogeneous 
spheres is a = b = 1/2, and that the restitution coefficient |(3]) may be justified, for 
vanishing speeds of impact, by computing the energy dissipated by ([99|l in a single 
compression-expansion cycle (see [3] chap. 3 and references therein). Although 
some of the claims in the literature appear to be questionable in their generality, 
because the vector field defining (j99| is non-Lipschitz for £ = and < b < 1, the 
problem is well-posed close to the position of static equilibrium, which is 

(100) i> = i> = i = 0, f = - ( 1 

That corresponds to the minimum of the energy 

(101) E = ^ + £. + -fL-\t\ a+a + T i> 

2 2 a + 2 
which is a non-increasing function of time: 

If we take a contact configuration (0, Xq, Yo, Yq) close enough to the static equilib- 
rium (|100p the system remains close to the minimum of the energy at all later times. 
As a consequence, we are guaranteed that X n and X n are in a neighborhood of 
—27 and 2/ry(7/p) 6 /( a+1 ), respectively. Then we follow the steps (08]) through ((5TJ) , 
using (|57|) as the time of flight, and we find that the asymptotic dynamics follows 
the map PI with a n -> p(37)~ 1 (7/p) b/(a+1) for 11 — > 00. Then we conclude that a 
sequence of repeated impacts, because of the deformability of the body, follows the 
same asymptotic law as the linear case. As a consequence, the restitution coefficient 
([3]) may be accurate only in the case of well-separated impacts, that is, when the 
time of flight between consecutive impacts is large enough to allow for dissipation 
of internal vibrations. 





9. Appendix 

Lemma 9.1. Given the map 

x n -\-\ — x n oi n x n -\- o(x n ) 
with < xq -c 1 and linin^oo a n = a > 0, then 

lim x n — and lim nx n = —. 

n — >oo n — >oo (X 

Proof. Since (a n — ot)x 2 n — o(x^), we may write 

(102) x n+ i = x n - ax n + o{x 2 n ). 

For small xo, the sequence generated by l|102p is positive, decreasing and IvaXn-,^ x n = 
0. For a given e, let us take N such that for any n > N 

(103) 1 - (a + e)x n < ^±1 < 1 _ ( a _ e ) Xn 

27 



We may assume xn < l/( a — e ) i then xjv+fc < l/(a — e); we prove that XN+k < 
l/k(a — e) implies XN + k+i < l/(k+ l)(a — e) and apply induction. Multiplying 
the right inequality in (|103p by 1 + (a — e)x n , we have 

^±I(l + (a- e )z„)<l 

Then 

, x N+lc 

XN+k+l < ~j j 7 \ 

1 + (a - e)x N+k 

where x/1 + (a — e)x is an increasing function of x and x^v+fc < l/k(a — e), then 



fc(ct-e) 



By induction 



then 



i.e. 



x n < 7 tttt r Vn > AT 



limsupriXn < Ve > 



By a similar argument we get 



lim sup nx n < — . 

n — >oo CV. 



lim inf nx n > — 

n — >oo q; 



which proves the proposition. □ 
Remark 9.2. It follows that for a map x n +i — x n — ax@ + o(x@) , with (3 > 1, 
(104) lim x n = and lim n^- 1 x n > 0. 

n — >oo n — >oo 

This can be seen by using the transformation z n — x^ 1 . We have 

= (1 - a^- 1 + o^- 1 ))^ 1 = 1 - aCS - ljxg" 1 + oixt 1 ) 

Zn 

We then apply Lemma l9~T1 to 

Zn+1 = Zn ~ a(/3 - l)z 2 n + o(zl) 

which gives Q104)) . 

Lemma 9.3. If x n > and x n +\ — x n + 0{x\) then J2 n x n = +°°- 
Proof. We have x„+i > x n — cx 2 n for a suitable c > 0. Hence, if we define 

2 

Zn-\-l — --■n CZ n 

z = x 

we have x n > z n for all n > 0. Applying Lemma l|9.ip to this map, the thesis 
follows. □ 

Lemma 9.4. If there exists a sequence {an}„ e pj with a n € (0, 1) that satisfies (E 
then lim T1 _ ) . 00 a„ = 0. 
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Proof. Let us define the functions / and g as 
(105) f(x) = 6x-Ax 2 + x 3 , 



, . , , 6x-8x 2 + 3x 3 

(106) g(x) = — -3 — . 

(1-x) 

The map i|6?j) is then written as 

(107) /(a„+i) =g(a n ) + 

with b n —>■ 0. Since / is bounded in [0, 1] and g tends to oo as x — ► 1, it is easy to 
see that there exists 5 > such that a„ < 1 — <5. We observe, that f'{x) > 0. With 
some algebra, we also find, for x € [0, 1), 

(108) /(z) - $(a;) < -- 



(1-x) 3 

In particular, the map ([67| has no fixed point in (0, 1). Let us assume that {a n } 
does not converge to zero. This means a n > e for infinitely many values of n and 
some e > 0. Inequality (|108j) implies that g(x) — f(x) < —e 2 for x > e. Let us fix 
m such that \b n \ < e 2 /2 for n > m and take any a n > e with n > m. Then the 
inequality a„+i < a„ would give, since / is increasing, 

b n = f(a n +i) - g(ctn) < f{a n ) - g(a n ) < -e 2 

which contradicts the choice of b n . Therefore a n+i > a n > e and, iterating this 
argument, the sequence a is eventually increasing and convergent to a number 
< a < 1 satisfying f(a) = g(a). Since this contradicts l|108p . the proof is 
complete. □ 
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